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Wo report results from searches for neutral Higgs bosons produced in pp collisions recorded by 
the DO experiment at the Fermilab Tevatron Collider. We study the production of inclusive neutral 
Higgs boson in the rr final state and in association with a b quark in the 6rr and bbb final states. 
These results are combined to improve the sensitivity to the production of neutral Higgs bosons 
in the context of the minimal supersymmetric standard model (MSSM). The data are found to 
be consistent with expectation from background processes. Upper limits on MSSM Higgs boson 
production are set for Higgs boson masses ranging from 90 to 300 GeV. We exclude tan /3 > 20 — 30 
for Higgs boson masses below 180 GeV. These are the most stringent constraints on MSSM Higgs 
boson production in pp collisions. 

PACS numbers: 14.80.Da, 12.60.Fr, 12.60.Jv, 13.85.Rm 



INTRODUCTION 



In the minimal supersymmetric standard model 
(MSSM) [1], the SU(2) symmetry is broken via two Higgs 
doublets; the first doublet couples to down-type fermions 
only while the second couples to up-type fermions. This 
leads to five physical Higgs bosons: two neutral CP-even 



*with visitors from "Augustana College, Sioux Falls, SD, USA, 
''The University of Liverpool, Liverpool, UK, ^UPIITA-IPN, Mex- 
ico City, Mexico, ''SLAC, Menlo Park, CA, USA, "^University 
College London, London, UK, ■'Ccntro dc Invcstigacion en Com- 
putacion - IPN, Mexico City, Mexico, SECFM, Universidad Au- 
tonoma de Sinaloa, Culiacan, Mexico, and ''Universitat Bern, Bern, 
Switzerland. 'DESY, Hamburg, Germany, tDeceased. 



bosons, h and H, one neutral CP-odd boson A, and two 
charged bosons . The neutral Higgs bosons are collec- 
tively denoted as (/). At leading order the mass spectrum 
and the couplings of the Higgs bosons are determined by 
only two parameters, conventionally chosen to be tan /3, 
the ratio of the two Higgs doublet vacuum expectation 
values, and Ma, the mass of the pseudoscalar Higgs bo- 
son. Radiative corrections introduce additional depen- 
dencies on other model parameters. Although tan /? is a 
free parameter in the MSSM, some indications suggest it 
should be large (tan^ > 20). A value of tan^ « 35 [2] 
would naturally explain the top to bottom quark mass 
ratio. The observed density of dark matter also points 
towards high tan /3 values [3] . 

At large tan /3, one of the CP-even Higgs bosons (/?, 
or H) is approximately degenerate in mass with the 
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A boson. In addition, they have similar couphngs to 
fermions, which are enhanced (suppressed) by tan/3 com- 
pared to the standard model (SM) for down-type (up- 
type) fermions. This enhancement has several conse- 
quences. First, the main decay modes become 4> ^ bb 
and 4> ^ T'^T~ with respective branching ratios S(<?!> 
bb) « 90 % and B{(l) r+r") w 10 %. Secondly, the 
main production processes at a hadron collider involve b 
quarks originating from the sea. Inclusive Higgs boson 
production is dominated by gluon fusion {gg(j)) and bb 
annihilation {bb(p), as shown in Fig. 1. The latter process 
may produce a b quark in the acceptance of the detector 
in addition to the Higgs boson. This associated produc- 
tion gb — > (j)b {bgbff)) is shown in Fig. Ic. In this case, the 
detection of the associated b quark is a powerful experi- 
mental handle for reducing backgrounds. 

MSSM Higgs boson masses below 93 GeV have been 
excluded by experiments at the CERN e+e~ Collider 
(LEP) [4]. The CDF and DO Collaborations have 
searched for MSSM neutral Higgs bosons decaying to tau 
pairs both inclusively [5, 6] and in association with a b 
quark [7]. The DO Collaboration has also searched for 
&(/) bbb production [8], which is challenging due to the 
high rate of multijet (MJ) production. Since these results 
have comparable sensitivities, combining them further 
enhances the potential reach. Recently, similar searches 
were performed at the LHC [9, 10]. In this Letter, we 
present a combination of three searches performed by the 
DO collaboration in the — )■ rr, bcp brr, and b(j) bbb 
final states. Since the inclusive and bgb<p production sig- 
nal samples in the di-tau final states are not mutually 
exclusive, the DO result presented in [6] can not be di- 
rectly combined with [7]. Hence, we re-analyse here the 
inclusive (jj ^ tt production: we require that there are 
no b jets, we extend the dataset to 7.3 fb^^ of integrated 
luminosity, and we increase the trigger acceptance and 
refine the treatment of systematic uncertainties. The di- 
tau channels are restricted to final states where one r 
lepton (r^) decays via r l-ii'fii'T and the other (r/,) 
decays hadronically. 

DETECTOR AND OBJECT RECONSTRUCTION 

The data analysed in the different studies presented 

here have been recorded by the DO detector [11]. It 
has a central-tracking system, consisting of a silicon mi- 
crostrip tracker and a central fiber tracker, both located 
within a 2 T superconducting solcnoidal magnet, with 
designs optimised for tracking and vertexing at pseudo- 
rapidities [12] \r]\ < 3 and [ryj < 2.5, respectively. A 
liquid-argon and uranium calorimeter has a central sec- 
tion covering pseudorapidities \r]\ up to w 1.1, and two 
end calorimeters that extend coverage to \r]\ » 4.2, with 
all three housed in separate cryostats [13]. An outer 
muon system, at \r]\ < 2, consists of a layer of track- 



ing detectors and scintillation trigger counters in front 
of 1.8 T toroids, followed by two similar layers after the 
toroids. 



TABLE I: Searches combined in this Letter. 



Final state 


£ (fb ^) Reference 


->■ Tfj,Th (&-jet veto) 


7.3 


b(j> hr^Th 


7.3 [7] 


b4> bbb 


5.2 [8] 



The integrated luminosities (£) [14] associated with 
each search are summarized in Table I. Di-tau events 
were recorded using a mixture of single high-pT muon, 
jet, tau, muon plus jet, and muon plus tau triggers. The 
efficiency of this inclusive trigger condition is measured 
in a Z — > r^T/i data sample with respect to single muon 
triggers. We also verify this measurement in a sample 
of TpT;, )+jcts events. Depending on the kinematics 
and on the decay topology of the Th, the trigger efficiency 
ranges from 80% to 95%. For the bbb analysis, we employ 
triggers selecting events with at least three jets. Most 
of the bbb data sample was recorded with ^-tagging re- 
quirements at the trigger level. The trigger efficiency for 
= 150 GeV is approximately 60% for events passing 
the analysis requirements. 

Muons are reconstructed from track segments in the 
muon system. They are matched to tracks in the inner 
tracking system. The timing of associated hits in the 
scintillators must be consistent with the beam crossing 
to veto cosmic muons. 

Hadronic tau decays are characterised by narrow jets 
that arc reconstructed using a jet cone algorithm with a 
radius of 0.3 [15] in the calorimeter and by low track mul- 
tiplicity [16]. We split the candidates into three dif- 
ferent categories that approximately correspond to one- 
prong T decays with no tt^ meson {th type 1), one-prong 
decay with tt^ mesons (r^ type 2) , and multi-prong decay 
{th type 3). In addition, a neural- network-based Th iden- 
tification (NNt) has been trained to discriminate light 
parton jets (u, d, s quarks or gluon) from hadronic r de- 
cays [16]. We select Th candidates requiring NNr > 0.9 
(0.95 for Th type 3). This condition has an efficiency of 
approximately 65% while rejecting ~99% of quark/gluon 
jets. 

Jets are reconstructed from energy deposits in the 

calorimeter [17] using the midpoint cone algorithm [15] 
with a radius of 0.5. All jets are required to have at 
least two reconstructed tracks originating from the pp 
interaction vertex matched within Ai?(track, jet-axis) = 
y/ {A-qY + (A(p)2 < 0.5 (where the azimuthal angle). 
To identify jets originating from b quark decay, a neural 
network 6-tagging algorithm [NNb] [18] has been devel- 
oped. It uses lifetime-based information involving the 
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FIG. 1: Main Higgs boson production mechanisms in the MSSM in the 5-flavor scheme where c and 6 quarks are included in 
parton density functions. The gluon fusion (a) and bb annihilation (b) processes dominate the inclusive production, while (c) 
is the dominant process for associated bcj) production. 



track impact parameters and secondary vertices as in- 
puts. 

The presence of neutrinos is inferred from the miss- 
ing transverse energy, fx- which is reconstructed as the 
negative of the vector sum of the transverse energy of 
calorimeter cells with |?7| < 3.2, corrected for the energy 
scales of all reconstructed objects and for muons. 



SIGNAL AND BACKGROUND MONTE CARLO 
SIMULATION 

Signal samples are generated with the LO event gener- 
ator PYTHIA [19]. The inclusive production is simulated 
with the SM ggcj) process. We checked that the kine- 
matic differences between bbcp and ggcj) do not have any 
impact on our final result. The associated production 
with a 6-quark is generated with the SM gb — ^ ipb pro- 
cess. The contributions to the b(j) cross section and event 
kinematics from next-to-leading order (NLO) diagrams 
are taken into account by using mcfm [20] to calculate 
correction factors for the pythia generator as a function 
of the leading b quark pt and r] in the range > 12 GeV 
and \r]''\ < 5. 

In the final states with a tau pair, the dominant back- 
grounds are due to Z — > TT(4-jets), diboson (PFW^, WZ 
and ZZ), H^-|-jets, ti pair and MJ production, the latter 
being estimated from data. Diboson events are simulated 
with PYTHIA while the Z+jets, W^-l-jets, and ti samples 
are generated using ALPGEN [21]. In the bbb channel, the 
dominant background is due to MJ production. We sim- 
ulate MJ background events from the bbj , bbjj, ccj, ccjj, 
bbcc, and bbbb processes, where j denotes a light parton, 
with the ALPGEN event generator. The small contribu- 
tion from ti production to the background is also sim- 
ulated with ALPGEN. The contribution from other pro- 
cesses, such as Z + bb and single top quark production, 
is negligible. 

The ALPGEN samples are processed through pythia 



for showering and hadronization. tauola [22] is used 
to decay t leptons and evtgen [23] to model b hadron 
decays. All samples are further processed through a de- 
tailed GEANT [24]-based simulation of the DO detector. 
The output is then combined with data events recorded 
during random beam crossings to model the effects of 
detector noise and pile-up energy from multiple inter- 
actions and different beam crossings. Finally, the same 
reconstruction algorithms as for data are applied to the 
simulated events. Data control samples are used to cor- 
rect the simulation for object identification efficiencies, 
energy scales and resolutions, trigger efficiencies, and the 
longitudinal pp vertex distribution. Signal, fl pair, and 
diboson yields are normalised to the product of their ac- 
ceptance and detector efficiency (both determined from 
the simulation) , their corresponding theoretical cross sec- 
tion and the luminosity. 



In the bbb final state, the relative contribution of the 
different MJ backgrounds is determined from data; its 
overall normalisation is constrained by a fit done in the 
final limit-setting procedure which exploits the dijet-mass 
shape differences between signal and background. In the 
di-tau channels, a dedicated treatment of the dominant 
Z ^ TT background has been developed to reduce its 
systematic uncertainties. The simulation of the Z boson 
kinematics is corrected by comparing a large sample of 
Z — > events in data and in the simulation. We mea- 
sure correction factors in each jet multiplicity bin as a 
function of the $* quantity introduced in Ref. [25], lead- 
ing jet r], and leading 6-tagged jet NNi,. This affects 
both the normalisation and the kinematic distributions. 
For the VF+jets background, the muon predominantly 
arises from the W boson decay while the Th candidate is 
a misreconstructed jet. The T4^-|-jets simulation is nor- 
malised to data, for each jet multiplicity bin, using a 
W{-^ /Lti/)-|-jets data control sample. 
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FIG. 2: Distribution of Mhat in the inclusive tt sample on (a) linear and (b) logarithmic scale, (c) T>f in the 6rr sample trained 
for = 100 GeV, and (d) for = 190 GeV adding the final requirements on Dmj and Dtt- AH th types are combined. The 
predicted signal is shown in the case of the m™'"' scenario {n = +200 GeV and tan/3 = 40). 



ANALYSIS STRATEGY 

In this section, we describe the search strategy as well 
as the selection of the final signal samples. Further details 
of the bTT and bbb analyses can be found in Refs. [7] 
and [8], respectively. 

Di-tau final states 

The TT and 6tt searches follow a similar strategy. 
We first define a common selection by retaining events 
with one reconstructed pp interaction vertex with at least 
three tracks, exactly one isolated muon and exactly one 
reconstructed Th- We require the muon to have a trans- 
verse momentum pt^ > 15 GeV, |ry''| < 1.6, and to be 
isolated in the calorimeter and in the central tracking sys- 
tem, i.e., Ai?(/x,jet) > 0.5 relative to any reconstructed 
jet. The Th candidate must have a transverse momentum, 
as measured in the calorimeter with appropriate energy 
corrections, p^ > 10 GeV, \r]rj < 2.0, AR{Th,n) > 0.5 
relative to any muon, and Th tracks must not be shared 
with any reconstructed muons in the event. The sum 
of the transverse momenta, p^^, of all tracks associated 
with the Th candidate must satisfy p^^ > 7/5/10 GeV, 
respectively, for Th types 1/2/3. We require the distance 



TABLE II: Expected background yield, observed data yield, 
and expected signal yields for the di-tau selections with their 
total systematic uncertainties. The signal yields are given for 
the m^=''= scenario {fi = +200 GeV and tan/3 = 40). 







TT 


bTT 


^(+jets) 




11547 ± 634 


218 ± 17 


tt 




25 ± 4 


183 ± 32 


MJ 




1343 ± 236 


36+6 


Other 




560 ± 25 


40+2 


Total back 


ground 


13474 ± 684 


476 + 40 


Data 




13344 


488 


Signal m0 


= 100 GeV 


1165 


81 


Signal 7710 


= 190 GeV 


70 


12 



along the beam axis between the Th and the muon, at 
their point of closest approach to the pp interaction ver- 
tex, Az{Th, A*) < 2 cm. In addition, the Th and the muon 
must have an opposite electric charge (OS) and a trans- 
verse mass Mt(m, I^t) < 60 GeV (100 GeV for type 2 ) 
where Mt{ijl,$t) = j2-pi^- $t- [l - cos A(^(/i, . 



7 



D0 , 5.2 f b"^ 
a) 3 jet 

Low-mass likelihood 



D0 Data 
Background 
Heavy flavor 




50 100 150 200 250 300 350 400 



M^^[GeV] 



D0 , 4.2 fb"^ 
b) 3 jet 

Low-mass likelihood 

—I— D0 Data 

Background 

Heavy flavor 
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FIG. 3: Distributions of the dijet invariant mass, taken from 
Ref. [8], in the signal region defined by Vtbt > 0.65 (a) and in 
a control region defined by Dbtt < 0.12 (b) for the dominant 
bbb 3-jet channel is shown. The line shows the background 
model, the solid histogram the component coming from bbb, 
the points with error bars show the data. The background 
is normalized to the data yield for illustration purposes. The 
difference between data and the background model is shown 
at the bottom of each panel. 



Inclusive tt selection 

For the inclusive tt selection, we tighten the require- 
ments on the Tfi transverse momentum to suppress the 



MJ background: p^'' > 12.5 GeV (15 GeV for th type 3) 
andpJi^'" > 12.5/7/15 GeV respectively for t,, type 1/2/3. 
We further reduce the W+jets background by requiring 
Mt(^, ^^t) < 40 GeV. We define Mhat, which represents 
the minimum center-of-mass energy consistent with the 
decay of a di-tau resonance, by 



Ml 



hat 



P. 



+ - 



Pt 



where E'^'^'^ is the energy of the /iT/j system and p^'^'' 
is its momentum component along the beam axis. We 
require Mhat > 40 GeV to suppress the MJ background. 
Finally, to prevent any overlap with the brr sample, we 
select only events for which no jet has NNi, > 0.25. 



brr selection 

The complementary sample with at least one 5-tagged 
jet with NNb > 0.25 constitutes the brr sample. This 
6-tagged sample suffers from large Z+jets, ti and MJ 
backgrounds. We build separate multivariate discrimi- 
nants, 2?MJ and V^j, to discriminate against the MJ and 
ti processes. We require Pmj > 0.1 and V^f > 0.1, then 
we combine NNb, I?mj, and P^j, to form a set of final 
discriminating variables Vf (one for each Th type and 
m^) to be used in the limit-setting procedure. Further 
details can be found in Ref. [7]. 



MJ background estimation 

In both di-tau channels, the MJ background is esti- 
mated from data control samples applying two different 
methods. The first is based on the small correlation be- 
tween the electric charge of muon and Th in MJ events. 
For each analysis, we select a data sample with identi- 
cal criteria as the signal sample but with the two Icptons 
having the same electric charge (SS). We subtract the 
residual contribution from other SM backgrounds from 
this MJ-dominated SS sample. We measure the ratio 
of the number of OS to SS events to be 1.09 ± 0.01 
and 1.07 ± 0.01, respectively, in the tt and brr chan- 
nels. We then multiply the SS sample yields by this ra- 
tio. This method is used in the inclusive tt channel but 
it suffers from large statistical uncertainties of the brr 
SS sample. Therefore, we develop an alternate method 
that uses a MJ-enriched control sample with identical 
requirements as applied to the signal samples but revers- 
ing the muon isolation criteria. In a MJ-dominated SS 
sample, obtained without any requirement on the num- 
ber of jets (A'jots), the ratio of the probabilities for a 
muon of a MJ-event to appear isolated or not isolated. 



-^iso/iiS = ^(AiisolMJ)/7'(^i^|MJ), is measured as func- 
tion of 77"^'' , pij!" , and leading-jet pr (if ^jots > 0). The ra- 
is then applied to the distributions of the non- 



tiO R; 



ISO / iso 

isolated- muon sample, predicting the MJ background in 
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the two signal samples. This method is used in the 
bTT study. In each analysis, the alternate method is 
used to determine the systematic uncertainty on the MJ- 
background normalisation. 

The distributions of Mhat for the tt study and two 
different Vf discriminants for the brr analysis are pre- 
sented in Fig. 2. The observed data, expected signal and 
background yields are given in Table II for the two di-tau 
event selections. 



bbb final state 

In the bbb analysis, at least three jets, each satisfying 
Pt > 15 GeV, |ry| < 2.5 and NNb > 0.775, are required. 
The two leading jets must have pr > 25 GeV. To im- 
prove the signal sensitivity, the events are separated into 
two channels, containing exactly 3 or 4 jets. The data 
and signal yields are given in Table III. In addition, 
a likelihood discriminant, I'bbb, based on six kinematic 
variables is employed. Two separate likelihoods, one for 
the mass region 90 < Ma < 140 GeV and the other for 
140 < Ma < 300 GeV, are used. The dominant heavy 
flavor multijet backgrounds are estimated using a data 
driven technique. The background in the triple b-tagged 
sample is estimated by applying a 2D-transformation in 
Mg^ and T>tbb, derived from the ratio of the number of MC 
events in the triple and double 5-tagged samples, to the 
double 6— tagged data sample. The method significantly 
reduces the sensitivity of the background model to the 
underlying kinematics of the simulated events and the 
modelling of the geometric acceptance of the detector. 
The appropriate composition of the simulated samples is 
determined by comparing the sum of the transverse mo- 
menta of the jets in each event in simulation and data for 
various ^-tagging criteria. The invariant mass distribu- 
tion of the jet pairing with the highest V^bb value is used 
as the final discriminant. The distribution for the domi- 
nant 3-jet channel is shown in Fig. 3a. In Fig. 3b, good 
agreement is observed between the data and background 
model in a control sample selected using an inverted like- 
lihood criterion Vbbb < 0.12. 



TABLE III: Observed data yield and expected signal yields in 
the bbb channel. The signal yields are given for the scenario 



Table II. 








bbb 


bbb 


iVjots 
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Data 
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166 
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FIG. 4: Comparison of the expected limits in the (tan /3, Ma) 
plane for the three channels separately, and their combination 
for the m™^'' scenario with (a) /i < and (b) /i > 0. 



Systematic uncertainties 

Depending on the source, we consider the effect of sys- 
tematic uncertainties on the normalization and/or on the 
shape of the differential distributions of the final discrim- 
inants. 

In the di-tau channels, the Z(-l-jets) background uncer- 
tainties are estimated using Z/j* — > IJ.~ data control 
samples, resulting in normalisation uncertainties of 3.2% 
(5%) for Z(-|-6-tagged jets) boson production, an inclu- 
sive trigger efficiency uncertainty of 3% (common to all 
simulated backgrounds) and a shape-dependent uncer- 
tainty of ^1% from the modeling of the Z boson kinemat- 
ics. The MJ-background uncertainty ranges from 10% to 
40% on the brr channel yields while it is found to be 
shape dependent in the tt channel (up to 100% at high 
Mhat)- For the remaining backgrounds and for signal, we 
consider uncertainties affecting the normalisation: lumi- 
nosity (6.1%), muon reconstruction efficiency (2.9%), Th 
reconstruction efficiency [(4-10)%], single muon trigger 
efficiency (1.3%), tt (11%) and diboson (7%) production 
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[GeV] 

FIG. 5: Constraints in the (tan/3, A/a) plane from the di-tau 
combination in the m™'"' scenario. These limits very weakly 
depends on the other MSSM parameters. 



cross sections. Further sources of uncertainty affecting 
the shape of the final discriminant are considered: the 
jet energy scale (10%) and the modeling of the ^-tagging 
efficiency (~4%) mostly affect the brr signal modelling 
but are negligible in the rr channel, while the energy 
scale (~10%) only impacts significantly the rr search for 
both Z boson background and signal Mhat distribution. 
With the exception of the Th reconstruction efficiency, Th 
energy scale and MJ estimation, which are evaluated for 
each Th type, these uncertainties are assumed to be 100% 
correlated across both di-tau channels. 

In the bbb channel, for the dominant MJ background, 
only systematic variations in the shape of the Mj^^ dis- 
tribution are considered, as only the shape, and not the 
normalisation, is used to distinguish signal from back- 
ground [8]. The dominant sources arise from the mea- 
surement of the rate at which light partons fake a heavy 
flavor jet and the ^-tagging efficiency. For the signal 
model, the ^-tagging efficiency (11-18%), the luminosity 
(6.1%) and the jet energy scale [(2-10)%] dominate the 
experimental uncertainties. 

Most of the experimental uncertainties are uncor- 
related between the di-tau and the bbb analyses with 
the exceptions of the 6-quark efficiency, luminosity, 
and jet energy scale, which are assumed to be 100% 
correlated. The theoretical uncertainties on the signal 
are other sources of correlated systematic uncertainty 
among all channels. They are dominated by parton 
density function uncertainties, renormalisation and 
factorisation scales. We assign an uncertainty of 15% on 
the theoretical cross sections that is correlated across all 
processes. 



RESULTS 

We combine the rr, brr and bbb channels using the 
modified frequentist approach [26]. The test statistic is 
a negative log-ratio of profiled likelihoods [27]: 

p(data|iio) 

where Hi is the test (background + signal) hypothesis, 
Hq is the null (background only) hypothesis and p are 
the profile likelihoods based on Poisson probabilities for 
obtaining the observed number of events under each hy- 
pothesis. We define CLg by CLg = CLs+b/CLb, where 
CLs+b and CL\j are the confidence levels for the test and 
null hypothesis respectively. We exclude signal yields 
with CLs < 0.05. 

The LLR quantity is computed from the M^at distri- 
bution for the TT channel, the Vf distributions for the 
brr channel and the M^j^ distribution for the bbb chan- 
nel. The NNLO SM cross sections dgg^ and dbb^ are 
taken from [28-35] and [36], respectively, while the NLO 
SM cross section abgb<p is taken from mcfm. The model- 
dependent MSSM to SM cross section ratios are com- 
puted with FEYNHIGGS [37]. To avoid double counting 
between the bbcf) and bgb4> processes, we obtain the ex- 
pected signal yield N^^^-^ in the di-tau channels by 

rr-f-X t X, ^ model i a ( ^ model ^ model\ 
^ - ^ ^ggi, + X [(Jbb<l> - ^bgb<t, ) 

\ A w ^ model 

+ ^bgb4> A (Tf,gf,0 , 

where the acceptances A are computed using the 
simulation and include the experimental efficiency. The 
two first terms of this equation refers to Higgs boson 
production without any b quark within the acceptance, 
while the third term is used for bgb(f> production. There 
is no difference in the experimental acceptance for the 
gg4> and bb4> processes with no outgoing b quark within 
the acceptance. Therefore, we set Abb4> = Agg^. The 
Higgs boson width, calculated with FEYNHIGGS, is also 
taken into account [8]. 

We test two MSSM benchmark scenarios [38], no- 
mixing and m^^^, and we vary the sign of the higgsino 
mass parameter, fi. The expected sensitivities for two 
^max scenarios are shown on Fig. 4 for the three different 
searches and for their combination. At low Ma, the brr 
channel dominates the sensitivity. For intermediate Ma, 
the TT and fcrr channels have similar sensitivities, while 
at high Ma, the bbb sensitivity becomes appreciable 
especially in /j, < scenarios. While the sensitivity in 
the TT+X channels are barely sensitive to other MSSM 
parameters than Ma and tan/3, the bbb signal yields is 
much more model dependent. Therefore we also provide 
a combination of the rr and 6rr searches only. We do 
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FIG. 6: Constraints in the {tan P, A^a) plane for difTerent MSSM scenarios from the combined Higgs bosons searches. 



not observe any significant excess in data above the 
expected background fluctuations and we proceed to set 
limits. The limit from the rr+X combination is shown 
in Fig. 5 and the full combination limits in different 
MSSM scenarios are shown in Fig. 6 . 



In summary, we present MSSM Higgs boson searches 
in three final states: rr, brr and bbb. These different 
searches are combined to set limits in the {tan (3, Ma) 
plane in four different MSSM scenarios. Furthermore, 
we combine the rr and brr channels to obtain MSSM- 
scenario independent limits. We exclude a substan- 
tial region of the MSSM parameter space, especially for 
Ma < 180 GeV where we exclude tan /3 > 20 - 30. These 
are the tightest constraints from the Tevatron on the pro- 
duction of neutral Higgs bosons in the MSSM and are 
comparable to the published LHC limits [9, 10], espe- 
cially at low Ma- 
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